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Abstract. We consider a primordial SU (2) vector multiplet during inflation in mod- 
els where quantum fluctuations of vector fields are involved in producing the curvature 
perturbation. Recently, a lot of attention has been paid to models populated by vector 
fields, given the interesting possibility of generating some level of statistical anisotropy 
in the cosmological perturbations. The scenario we propose is strongly motivated by 
the fact that, for non-Abelian gauge fields, self-interactions are responsible for gener- 
ating extra terms in the cosmological correlation functions, which are naturally ab- 
sent in the Abelian case. We compute these extra contributions to the bispectrum of 
the curvature perturbation, using the SN formula and the Schwinger-Keldysh formal- 
ism. The primordial violation of rotational invariance (due to the introduction of the 
SU{2) gauge multiplet) leaves its imprint on the correlation functions introducing, as 
expected, some degree of statistical anisotropy in our results. 

We calculate the non-Gaussianity parameter /jvl, proving that the new contributions 
derived from gauge bosons self- interactions can be important, and in some cases the 
dominat ones. We study the shape of the bispectrum and we find that it turns out to 
peak in the local configuration, with an amplitude that is modulated by the preferred 
directions that break statistical isotropy. 
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1. Introduction 



The standard cosmological model is based on tlie idea of inflation as driven by a single 
slowly-rolling scalar fleld whose primordial quantum fluctuations are responsible for both 
CMB perturbations and the large scale structures of the universe [1]. Current observa- 
tions indicate that the universe on large scales is homogeneous and isotropic. The degree 
of deviation from this smooth background is provided by the temperature fluctuations of 
the CMB which are of the order of 10~^ and are almost scale-independent [2, 3, 4, 5, 6]. 
The standard inflationary scenario predicts the observed power spectrum of cosmic fluc- 
tuations. On the other hand, there is still a lot of room for alternative scenarios without 
ruining the current agreement with observations. Alternatives include, for example, 
multifleld models [7, 8, 9, 10, 11, 12, 13], the curvaton scenario [14, 15, 16, 17, 18, 19] 
and theories with non-canonical Lagrangians as k-inflation [20, 21], DEI inflation [22, 23] 
or ghost inflation [24]. So far, observations have not provided us with the precision mea- 
surements which would allow to discriminate among all these different models, but new 
projects, such as the Planck [25] satellite which was just launched, can potentially reach 
the levels of precision necessary for this purpose. One important feature of the CMB 
anisotropies to be decoded is the degree of non-Gaussianity [26, 27]. A random fleld is 
deflned as Gaussian if all the information is contained in its two-point correlation func- 
tion. In single-fleld, slow-roll inflation the level of non-Gaussianity is very small since 
the connected parts of higher order correlation functions are proportional to powers of 
the slow-roll parameters [28, 29, 30]. If a larger non-Gaussianity will be detected in 
the future, this could then rule out the minimal model. Recently another alternative to 
the standard inflationary scenario has attracted some attention where primordial vector 
flelds play a non-negligible role during inflation [31, 32, 33, 34, 35, 36, 37, 38, 39, 40] . 
An interesting possibility is that these extra flelds might be responsible for generating 
an observable level of non-Gaussianity [36, 38, 35]. This is not the only motivation to 
consider such models. First of all vector flelds are present in particle physics models and 
that suggests that one has to include them in a theory of the early universe. This can 
be done without necessarily affecting current observations, for example keeping their 
contribution to the energy density of the universe way below the one from the infla- 
ton. However, a peculiarity of the vector flelds, which the scalar flelds do not possess, 
is to violate rotational invariance as some observations seem to suggest. Violation of 
rotational invariance would affect the correlation functions introducing some degree of 
statistical anisotropy. The power spectrum is constrained in this sense if a single pre- 
ferred direction of space is involved [41]. This is something to keep in mind for example 
when constucting models of vector inflation. The isotropic background expansion can 
be preserved to a good approximation in different ways, such as by keeping the energy 
density of the vector subdominant w.r.t. the one of the inflaton, or by considering the 
existence of many randomly oriented vectors [33] . Nevertheless, the presence of the vec- 
tor flelds introduces statistical anisotropy in the power spectrum and in the bispectrum 
as well [36]. On the other hand, no observational constraints are currently available for 
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the two-point function if more than one preferred direction is involved. Similarly, for 
higher order correlators, no obervational constraints of this kind are available. 

The effects of a primordial vector field on the power spectrum and bispectrum of 
the curvature fluctuations ( were recently investigated in [36, 38] using the 6N formal- 
ism [42, 43, 44, 45] both in vector inflation and in the vector curvaton models, and also 
in [35] for a model of hybrid inflation. 

Up to now only primordial Abelian vector fields were considered. In this paper we 
wish to explore a more realistic model of gauge interactions of an SU (2) vector multi- 
plet. The reason for considering a non- Abelian theory relies on the fact that, differently 
from what happens in the Abelian case, the vector field components are involved in 
self-interactions as described by cubic and quartic terms in the Lagrangian. Gauge 
self-interactions naturally produce additional contributions to the bispectrum w.r.t. the 
Abelian case. The computation of these new contributions represents the main purpose 
of our work. 

The paper is organized as follows: in Section 2 we introduce the SN formula for 
the curvature perturbation and the Lagrangian of our system; in Section 3 we carry out 
the calculations of all the terms that do not require the intervention of the Schwinger- 
Keldysh formula; in Section 4 we derive the contribution to the bispectrum that arises 
from gauge bosons self-interactions; in Section 5 we study the shape of the bispectrum; 
in Section 6 we estimate the non-Gaussianity parameter f^L, in Section 7 we draw our 
conclusions. We include four Appendices: A, on background and first order perturbation 
equations for the gauge multiplet; B, where we derive the expression of the number of 
e-foldings of inflation and its derivatives in the presence of vector fields; C, which collects 
lengthy expressions for functions appearing in the final results for the bispectrum; D, 
where we give some details about the dependence of the bispectrum from the angles 
between wave and gauge vectors in a sample spatial configuration. 

2. Bispectrum of the curvature perturbation in the 6N formalism 

In the the SN formalism the curvature perturbation ({x) at a given time t can be 
interpreted as a geometrical quantity indicating the fluctuations in the local expansion 
of the universe; in fact, if N{x,t*,t) is the number of e-foldings of expansion evaluated 
between times t* and t, where the initial hypersurface is chosen to be flat and the final 
one is uniform density, we have 

C(f, t) = N{x, t*,t)- N{t\t) = 5N{x, t). (1) 

The number of e-foldings N{x, t*,t) depends on all the fields and their perturbations on 
the initial slice. In principle, since the fields are governed by second order differential 



3 



equations, it should also depend on their first time derivative, but if we assume that slow- 
roll conditions apply, then the time derivatives will not count as independent quantities. 
In a theory which includes vector fields, spatial isotropy is not expected to be preserved 
and therefore the power spectrum and the bispectrum can in principle depend both on 
the length adn direction of the wavevectors ki, rather than just on their moduli 



(2) 
(3) 



In the next section, we are going to spell out the 6N expansion in terms of the fields of 
our theory; later we will replace it into Eq. (3) and calculate the bispectrum to tree-level. 



2.1. Specializing to our theory 

Our theory includes a scalar field (p playing the role of the infiaton and an SU{2) gauge 
multiplet (a=l,2,3) non-minimally coupled to gravity. The metric of an FRW flat 
spacetime ds"^ = dt^ — o?5ijdx^dx^ is employed and held unperturbed. The action is 



S 



d X\ 



m%R 



1 /lf2 

^ a=l,2,3 ^ a=l,2,3 



(4) 



where La is the scalar field Lagrangian and 



FL = d,B^,-d,B; + g,e^''^B';^B: 



(5) 



(we name gc the gauge coupling). Notice that B'^ are the comoving fields, the physical 
fields being given by = {Bq, (l/a)-B"). The coupling between the gauge bosons and 
gravity is "hidden" in the effective mass M: by definition = ttlq + ^R, where mg is 
the gauge bosons mass and ^ is a numerical factor; using the slow-roll approximation, 

2n 



the Ricci scalar is i? = — 6 
The quantum fluctuations 



7 + 



-12H' 



or the vector fields can be expanded in terms of creation and 
annihilation operators (this is always correct if we are dealing with free gauge boson, 
their background equations of motion being of Klein-Gordon type) 

d^q 



where 



(2vr) 



a, A 



X=L.RI 



a ,A 



(6) 



(7) 



dr] = dt/a{t) is the conformal time and A labels left, right and longitudinal polarization 
states. 

The perturbation of the number of e-folds on large scales is determined by the 
perturbations of the fields at the initial time and the complete 6N formula then reads 
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C(f , t) = N^6<p + N^SA';^ + -N^^ {5ct>f + -iVjM^M^ + N^M^A^ + ... (8) 
where 



--(f).. ^Hi^). 

and so on. is the number of e-foldings between the time t* on the initial hypersurface 
and the final time rj at which the bispectrum of the C is being evaluated. A convenient 
choice, that we will follow for the rest of the paper, is to take rj* to be close to the 
time when a given perturbation mode of wavenumber k crosses out the horizon. The 
coefficients in Eq.(9) will be calculated in Appendix B. Notice that the e-folding numbers 
and their derivatives are defined w.r.t. the physical fields, whereas the action (4) is 
written in terms of the co moving fields B^. 

Plugging the 6N expansion into Eq.(3), we get an infinite series of terms. Retaining the 
tree-level ones (i.e. the terms that are formally equivalent to the product of two power 
spectra), we are left with 



In the previous expansion we have not included terms that would require some kind of 
coupling between scalar and vector fields, i.e. terms that are proportional to N^N^Nl, 
N^Ni^, N^NiNi^, N^Ni, N^a^ic and N^N^^Nl Also, the terms that are proportional 
to A^'q = [dN) / [OAq) do not appear in the previous equation since the Aq fields are 
set equal to zero. The reason for this will be explained in details in Appendix A (see 
discussion after Eq (A. 7)); here we anticipate that, with a zero expectation value for the 
temporal components of the background gauge fields and with some upper bounds on 
the SU (2) coupling constant of the theory, it is possible to obtain a slow-roll evolution 
for the fields during inflation. 

It is useful to notice that the terms in the first four lines of (10) are also present in the 
Abelian case, while the one in the last line is strictly non-Abelian: it provides a non-zero 
contribution only thanks to gauge bosons self-interactions. 
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3. Power spectrum and the 'Abelian' terms of the bispectrum 



It is useful to write the expression for the power spectrum of the curvature perturbation 
in the presence of a vector muhiplet. Let us define 

P^^ = (l/2)(P]J^±Pf), (11) 

where 

P^' ^5at5A%ik,t*)5A'^{k,t*), (12) 

=5ab6Alik,t*)6A'*ikt*), (13) 

P^i = Sat6Al^,{k, e)5At^^{k, n (14) 

{6 A indicating the eigenfunctions for the gauge fields). The power spectrum of C can 
then be written as 



Pdk) = N^P^ik) + KNl [5,,Pf + hk, {P,t.^ - Pf 



i^ijkkkP- 



ab 



P''°{k) 1 + g"^ [k- Na){k- Nb) + is^'k ■ [Na x Ni 



^abr 



where 



^ab 



pab pab 

long + 



mPs+ N,-N,)pf 



s = 



N^P^ + [N, ■ Nd) Pf 
and we have factorized the isotropic part of the power spectrum 



(15) 

(16) 
(17) 



P''\k) = N^P^ik) + [n, ■ N,) Pf . (18) 

In Eq. 15 g"''' and s"^ represent the amplitude of statistical anisotropy, along the preferred 
directions specified by the vectors Na = and A^;,. Eq.(18) can also be written in the 
form 



P''\k) = NlP^ 
introducing the parameter 



cd' 



P^ 



(19) 



'cd 



(20) 



Let us now move to the bispectrum. The first three lines of Eq.(lO) can be easily 
evaluated 
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Bcih, k2, hf''^ D \nIN^^ [P5^{h)Ps^{k2) + perms. 



+ -KNIN':^ [YVlt{h)Y\'f,{k2)+ perms. 



1 

2' 



;N^KNl, [Ps^{k,)YVt^{k2)+ perms. , 
where we have defined 

n^^(^) ^ T:;'^-{k)Pi' + iT°f\k)Pi' + Tlp'ik)Pit^. 

In the previous equations, Tf.*"^", T"'^'^ and T/°"^ are defined as in [36] 



(21) 



(22) 



with e^{k) 
e\k) = k = 



T^ri^) = ^hk)ef{k) + ef (^)ef (^), 
(^) ^^\e^{k)ef{k) - ef{k)ef{k) 



T!p<^{k)^eme*/{k) 



(sin 9 cos 0, sin 6* sin 0, cos 6') , so 



(23) 
(24) 
(25) 



^(cos6' cos — i sin 0, cos6' sin + i COS0, — sin ^), e^{k) = e*^{k) and 



^ij (^) ~ ^"ij ~ kikj, 

'Pij ^ {k) = kikj. 



(26) 
(27) 
(28) 



Notice that if gauge indices are suppressed in Eq.(21), the Abehan expression is recovered 
(see [36]). 

It is convenient to rewrite the expression in the second line of Eq.(21) as follows 



where 

M\k) = Pl\k) [Na + P^^kYk {k ■ N,) + iq''\k)k X 
and we define 



pac TDac 

^long + 



PI 



pa 



PI 



Similarly, the third line in Eq.(21) becomes 



(29) 

(30) 

(31) 
(32) 



N^NlNl,P^{k{)Ilf;{k^) = N^NMk{)Mt{k^) 



(33) 
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The contribution from the fourth hne in Eq.(lO) can be calculated using the 
Schwinger-Keldysh formalism [46, 47, 48]. It is a well-known result tha the sum of 
this last term and the one from the first line of Eq.(lO) is proportional to the slow-roll 
parameters [28, 29]. 

It is important to notice that, in models where -P^^^g = -P"^ and P"^ = 0, the coefficients 
gab^ ^ab^ ^ab ^ab equal to zcro and therefore both the power spectrum and 

the Abelian part of the bispectrum of C, become isotropic. In general P"^ = 0, unless 
there is some mechanism producing parity violation for Ai. The condition P^^"^ = P^^ 
is more subtle. It could be realized if the longitudinal and the transverse parts of the 
vector fields evolved in the same way [36]. 



4. Calculation of the 'non- Abelian' terms of the bispectrum 

The last line of Eq.(lO) also requires the Schwinger-Keldysh formula 



{n\e{t)\n) 







T e 



ej{t) 



T e 



-iJlHi{t')dt' 







(34) 



where 6(t) is a field operator, \Q) represents the vacuum of the interaction theory, 
T and T are time-ordering and anti-time-ordering operators. All the fields are in the 
interaction picture, as the subscript / indicates. 

Our next step will then be to write explicitly the action for the gauge bosons and look 
at the detailed expression of their interaction Hamiltonian. Using the definition of F^^,, 
the action for the gauge fields becomes 



(35) 



- g,e^'^g^^g''P{d,B:)BlBl - ^gy^e'^'^g^^^g'^PBlBlBiBl - ^g^'^B^B^, 



The interaction Hamiltonian to third order in the gauge field perturbations is thus 
made up of two contributions coming respectively from the third and fourth order 
interations in the lagrangian 



Hint = gce''^''(d^SB^)SB^^'6B"' + gy^e'^^BlSBlSB^^'SB'^'' . (36) 
The correction to the bispectrum (10) due to these interactions has the form 



+perms. + c.c. 



(37) 



where 



Iit^{k) = 7;7^"(fc)Pf + iT^f\k)Pl'' + Tl°'''{k)Pi 



ab 
long 



(38) 
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and P± = (1/2)(Pr ± P^), being defined as the product of tlie eigenfunctions 
6BR{k,r]*) and 6B}^{k,ri), and so on for Pi and Piong- Notice tliat 77 is tlie time at 
which the three point function for ( is being evaluated. Also, we dropped the gauge 
indices in P, since the perturbations have the same time and momentum dependence 
for the different gauge fields. This last point can be discussed looking at Eq. (A. 15) for 
the perturbations of SB"", the equations for the free fields (which are the ones we need 
to compute the power spectrum and the bispectrum) are obtained from (A. 15) setting 
gc = 0, i.e. suppressing the interaction terms, and thus we get the same equation for all 
the components of the gauge multiplet. 



4.1. Bispectrum from third- order interactions 

Let us begin with the bispectrum contribution from Eq.(37) at the lowest order in the 
couphng gc 



X r dr{a\,i)IVt^tl'ljlt9''^g''^ikf^g, + perrns. + c.c m 
J —00 



where g^^ is the unperturbed background metric. The factor ik[^^ originates from 
the spatial derivative in the interaction term considered, whereas come from ^/—g 
together with having moved from cosmic to conformal time dt = adrj. In order to write 
Eq.(39), we used the third order interaction Hamiltonian 



H. 



(3) 
int 



_abc 



g'^g^' (diSB^) 6Bl6B^ + g^^ {d^SB^) 5B]5Bl + g'^ {do6B1) 6B',6B^] . (40) 



The calculation of the tree-level contributions due to the last two vertices in Eq.(40) 
will not be carried out since they contain the time components SBq, which turn out to 
be affected by instabilities. In fact, SBq is a non-dynamical mode that, for the free-field 
case, can be expressed in terms of the longitudinal component as follows [31] 



dt (k ■ 6B) 

The longitudinal modes of the theory have negative kinetic energy and are therefore 
considered unstable. In addition to this, their equations of motion suffer from 
singularities around the horizon crossing time. However, a regular solution was provided 
in [36], which we are going to employ when handling these modes (Eq. (58)). 
The factor l/a^{ri*) outside the integral accounts for the fact that the action is written 
in terms of the comoving fields B^ whereas the bispectrum of ( is written in terms of 
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the bispectrum of the physical fields A'}. 

In our model, -P_ = (i.e. Pr = Pl), therefore the integrand of Eq.(39) can be written 
as the sum of products of P+ = Pr and Piong weighted by products of coefficients T^^'^^ 
and T'°"3 



rpEEE 

ijk 



ki Sikkji — kiikjikki — ki^kjikk^ — 6ik{ki ■ ^2)^2 + (^1 ■ ^2) (^i3^j2^fc3 



+ kiikj2kki) + {ki ■ ks)kiikjikk3 - {ki ■ k2){ki ■ k3)kiikj2kk3 



nEEl 



rpElE 

ijk 



rplEE 

ijk 



- ki [kiskjikks - ki3kj2kk3{ki ■ /c2) - kukjikksiki ■ ks) 
+ kiikj2kk3{ki ■ k2){ki ■ h) 

= ki 5ikkj2{ki ■ h) - ki3kj2kk3{ki ■ h) - kiikj2kki{ki 
+ kiikj2kk3{ki ■ k3){ki ■ k2) 

= ki kiikjikki — kiikjikk3{ki ■ ks) — /cii^j2^fci(^i ■ ^2) 
+ kiikj2kk3{ki ■ k2){ki ■ h) 

= ki kiikj2kki{ki ■ h) - kiikj2kk3{ki ■ k2){ki ■ h) 
= ki \kiikjikk3{ki ■ h) - kiikj2kk3{ki ■ k2){ki ■ h) 
= ki [ki3kj2kk3{ki ■ h) - kiikj2kk3{ki ■ k2){ki ■ h) 
= ki kiikj2kk3{ki ■ k2){ki ■ h) 



(42) 

(43) 

(44) 

(45) 
(46) 
(47) 
(48) 
(49) 



rpUE 

ijk 

rplEl 

ijk 

rpEll 

ijk 

rplll 

where where we used the notation T^^^ = T[i"^"'Tj^"'Tl°^^6'^"^6''^ks and so on. 
The next step before carrying out the final integration over time in Eq.(39) is to sum 
over all the permutations separately for each one of the eight Tj"^''^ coefficients. Because 
of the antisymmetric properties of e"'^'^, the sums provide results of the form e"'^'^Sab, 
where Sab is a symmetric tensor. We thus find out that, independently of the results 
of time integrals (i.e. independently of the particular form of the wavefunctions for 
transverse and longitudinal modes), the contribution to the bispectrum from third- 
order interactions is zero. 

Let us then move to considering the quartic interaction. 



4.2. Bispectrum from fourth-order interactions 
The contribution from the quartic interaction is 



iCkSk^Ck,) ^ i^^m E ^^)^c4f^) K'KN'c'Ar"''^ 



i=l,2,3 H\k 



3 

k Ae^aa'c'ab'e 



X J2([d^) Tt^T^Xi + perm.. + c.c. 



(27r)353( ^ k,)gl^{Hx*fY.K.^i + P^^^s. (50) 

a/37 
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1=1,2,3 a/37 



with X = —krj. (Hx*)^ comes from the l/a^{r]*) factor outside the time integral 
and the 1/H factor comes from having a factor of a inside the time integral, i.e. 
/ drja^r]) = — J dx/{Hx). This factor originates from the background field appearing 
inside the integral (the interaction term is quartic and can be developed as a product 
of three perturbations times a background field), i.e. -Bf (77) = a{ri)A1{ri), where A is 
the slowly-rolling physical field (and hence taken out of the integral). Notice that, once 
again, we are setting Aq = 0, so terms proportional to do not appear in the formulas. 
The indices 7 as usual run over the longitudinal and the transverse components. 
(/ dx)^p^ is a shorthand for the integral over the wavefunctions associated with the 
internal legs times the wavefunctions associated with the external ones. The definition 
of the symbol K^ii^ can be read from the first two lines in Eq.(50). It is important to 
sum over all of the permutations, i.e. over all the possible field contractions, which we 
are going to list below 



cf ^ = TtpTf,T2p5^,,6yA'ci (51) 

cf ^ = TrpT^pT^MS,,,6,, (52) 

cf^ = T,^TfpTl6.,Mc'd (53) 

cf ^ = TtiT^pT2p5a'A'd5c'b (54) 

cf ^ = TtpT^pTiJ^Sa'Mc'c (55) 

cf ^ = TtpT^iTl6.,,6,,J,,,. (56) 

Let us write the expression of the wavefunctions of the gauge fields. We are going 
to work in the regime where = —2H'^, i.e. we take mo -C H and = 1/6. This 
is the situation where the transverse mode of the 6B field is governed by an equation 
of motion similar to the one of a light scalar field in unperturbed FRW background 
(see Eq.(A.16)). The solution is very well known and, after matching with the vacuum 
solution at early times, is given by 

-^3/2(3;) + iJ^3/2{x) . (57) 



2^/k 

where x = —krj. We have dropped the gauge index since the equation of motion is unique 
for all of the gauge bosons (see Eq.(A.16) and (A. 17) and discussion below Eq. (38)). 
As to the longitudinal mode, there is an on-going discussion about its instability [36, 
37, 39, 40]; according to some authors [37, 39], there are two sources of instability: the 
first one is due to the fact that the equation of motion (A. 17) is singular for a given 
mode k equal to the effective mass of the vector field; the second one arises from the field 
behaving like a 'ghost', i.e. having a negative energy. In spite of the divergence affecting 
the linearized equation of motion, the authors of [36] provided two independent exact 
solutions to (A. 17). The initial conditions were formulated after rescaling the vector 
field with the introduction of 5B = (a\M\/k)SB^^, which allows for the well-known 
initial conditions 6B = e~*^^/ ^/2k at very early times. We are going to use the solution 
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and the initial condition prescription of [36], so the longitudinal mode function acquires 
the following form 



2Vk \ X J 
The terms Kais-y below (permutations are included) then read 



(58) 



Keee 
Kin 

KuE 

Keei 



IEEE 



2Ak^klklklx*^ 



[Aeee + {Beee cos X* + Ceee sin x*) EiX* 
[Am + {Bm cos x* + Cm sin x*) EiX*] 



192k^kfkiklx*^ 

"^"^ — [AiiE + (BiiE cos X* + Cue sin x* ) EiX* 



9Qk^klklklx*^ 
Ieei 



[Aeei + {Be El cosx* + C eei sin x*) EiX* 



(59) 
(60) 
(61) 
(62) 



T ^aa'b' ^ac'e 

J-EEE = £ £ 



+ 



ASk'^klklklx*'^ 

where A, B, C and D are functions (to be provided in Appendix C) of x* and of the 
momenta ki = \ki\, i = 1, 2, 3 and E^ is the exponential-integral function. 
The factors Ia/3'y are given by 

6 (iV"' ■ iV"') (iV^' ■ A^) 

+ (iV''' ■ A^) [( - 2 (fcs ■ iV"') (A;3 ■ N'') - 2 (ki ■ N^') {ki ■ N"') 

+ [h ■ N"') {k ■ N'') k ■ k + (^3 ■ iV"') [h ■ N"') k ■ k) + (1 ^ 2) + (3 ^ 2) 

"(2 (iV'^' ■ N'') [k ■ N'') {h ■ A')) + (2 ^ 1) + (2 ^ 3)" 

(k ■ A'[2 {k ■ N"') {h ■ N'') {k ■ N'') +2[k- N"') {k ■ N"') {k ■ N'') 

- {k ■ N'^') {k ■ N'') [k ■ N^') k-k- {k ■ N'^') {k ■ N"') [k ■ N^') k ■ k: 

+ (2 ^ 1) + (3 ^ 2)]" 

= (^ki ■ N""') [k ■ N"') [k ■ N"') {k ■ k) {k ■ A') 

- [k ■ iV"') {k ■ N"') [k ■ iV^') (^1 ■ k) {k ■ A j ) + (1 ^ 3) + (2 ^ 3)" 

"(2 {k2 ■ A') {k ■ iV"') [k ■ N"') {k ■ N"') ) + (1 ^ 2)" 

( {{k ■ N"') {k ■ N"') + {k ■ N"') {k ■ N"')) {k ■ N'') {k ■ A') k ■ k 

+ (1 ^ 3) + (2 ^ 3)]" 

4 (iV^' ■ 1") {ks ■ N"') [k ■ N'') 

k ■ N"') [k ■ A!^) {{k ■ N"') {k ■ N'') + {k ■ N"') {k ■ N"')) k ■ k 

+ (2 ^ 1) + (2 ^ 3)" 

2 [k ■ A') {k ■ iV"') [k ■ N'') [k ■ N'')) + (1 ^ 2) + (2 ^ 3) + (1 ^ 3) 
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nil 



niE 



(63) 

(64) 



+ 



(65) 



T ^aa'b' ^ac'e 

J-EEl = ^ ^ 

+ 



+ 



(iV«' . iV''') ■ N"') [k ■ A^) + (iV"' ■ iV^') (^3 ■ iV"') (^k (66) 



where i ^ j means 'replace ki with kj\ whereas i ^ j means 'exchange ki with kj\ 
The terms that were not included in the hst above, i.e. Kiee, Keie, Keiu Kieu are 
derived from Keei and Kue by an appropriate exchange of the momentum labels, for 
example Kiei and Keie are obtained respectively from Kue and Keei by permuting /c2 
with /ca, Keu and Kiee from Kue and Keei by permuting ki with k^. 
As one can see through Eq. (59) to (66), the resulting bispectrum is made up of isotropic 
parts that depend only on the moduli of the wavevectors ki and that are modulated by 
the coefficients lap'y- These coefficients contain the information about the breaking 
of rotational invariance, depending on the angles between the wavevectors and the 
preferred directions A^". 

Given the fact that the instability of the longitudinal mode of vector perturbations is 
still considered a controversial issue, a possible alternative to the exact solution provided 
in [36] is to parametrize the longitudinal mode in terms of the transverse solution as 
follows 

5B\\ = n{x)6B^ (67) 

where n is an unknown function of x = —kr]. The most relevant contribution to 
the time integrals in Eq.(50) is due to the region around horizon crossing x* (see for 
example [29]), therefore a good approximation when integrating can be to treat n as a 
constant evaluating it at x = x*. The sum over K^is-y then becomes 

E = ^y^^ \Ieee + n^ix*) (Ieei + Ieie + Iiee) + n'^ix*) {Iue + Iiei + Ieu) + n^{x*)Iui 

a/37 

(68) 

The accuracy of this approximation has been tested considering again the solutions 
given in Eq.(57) and (58), which satisfy the relation = \/25B^ at horizon crossing. 
If we set n{x* ) = ^2 in Eq.(68) and compare each one of the terms in round brackets 
with the corresponding exact time integrals provided in Eq.(59) through (62), we get a 
factor of 2 to 5 between the exact time integrals and their approximations for the most 
part of the (x2,X3) range. 

The particular case n{x*) = 1, corresponding to logitudinal and transverse modes which 
evolve in the same way, is worthy of mention: the anisotropic part of Eq.(50) cancels 
out leaving an isotropic contribution only 



j=i,2,3 n \ n / 

Keee , 

X — h perms. (69) 

Ieee 

The way this result is achieved can be easily seen from the expression of the T°'T^T'^ 
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products, from which Kap^ were evaluated. We hst these products for one of the 
permutations 

^u^fpT^p = SiiSjp + 5ii {—kjskks — kj2kk2 + %2^fc3^23) — kii (Sjkkii — ky^^ — ki22 + ^123^23) 
^jp^fcp = ^ii {kj^kk-i - kj2kk3k23) - hi (/cisa - ki23k23) 

^uTjpT^p = Sii {kj2kk2 — /Cj2^fc3^23) ~ hi {ki22 — /^123^23) 
"^il^jp^kp — hi {^jkkil — A;i33 — ki22 + ki23k23) 

T^iTjpT^p =hi{ki22- A;i23^23) 
TiiTfpTlp = kii {ki3s - ki2'ik2'i) 
T^fT-pTlp = 6iikj2kk3k23 - A;/iA;i23^23 
TiiTjpTlp = hiki23k23 

where kabc = kiakj^kkc and kab = ka ■ k^ {a,b,c = 1,2,3 running over the external 
momenta). Notice that summing over the eight terms in the equations above only 
leaves SuSjp. 

5. Shape of the bispectrum 

As we saw in Eq.(50), the bispectrum in the presence of an SU{2) gauge multiplet can be 
written in terms of isotropic parts, i.e. functions of x* and of the moduli of the external 
wave vectors fcj, weighted by anisotropic coefficients /q,/3^ which depend on the angles be- 
tween the (wave and gauge) vectors. One way of studying the profile of the bispectrum 
is to analyse the isotropic parts separately in order to understand what is the preferred 
configuration, for example if it resembles the local or the equilateral shapes [52]. Once 
the answer to that is found, we calculate the anisotropic coefficients in that specific con- 
figuration. Let us employ the variables X2 = /c2/^i and X3 = k-s/ki as in [52]. Plotting 
the isotropic part, one can notice that for most of the functions KajS-y/Ia/B-y, which sum 
up to provide the bispectrum in (50), the regions with the highest values are around 
X2 = 1, X3 = (see plots in Fig.l). The graphs that don't have their maxima in this 
region, i.e. ruE and teie, which peak around X2 = X3 = 1, and tiee, which peaks around 
X2 = X3 = 0.5, show negligible amplitudes w.r.t. the previous graphs. Therefore we can 
safely state that the bispectrum peaks in the 'local' configuration k2 ~ ki, k^ <C ki. 
In fact this result can be understood recalling that the contribution to the bispectrum 
comes from the quartic interaction terms in Eq. (36), which are local in space, and do 
not involve gradient terms of the fields. 

We can now calculate the anisotropic coefficients in the local configuration. We have 
the freedom to pick up a coordinate system (x, y, z) with ki and k2 pointing along z and 
k^ along X. In the general case of three differently oriented vectors Na-, the coefficients 
/q,/3^ depend on a number of angular variables that is bigger than two and, therefore, 
they are not easy to plot. What we can do is looking at some particular cases with 
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Figure 1. Plot of Vap-y = &{x2 - a;3)9(a;3 - 1 + X2)x2xlRa[3'y{x2, X3) , 
where we define Rap-^ — k\{Kaf}j/ lapj)- Ttie Heaviside step functions 
Q help restricting the plot domain to the region {x2,X3) that is allowed 
for the triangle fci + ^2 + fca = (in particular, we set X3 < a:;2).We also 
set X* = 1. -,r 



specific orientations of tlie A^^ w.r.t one anotlier and to tlie coordinate frame we built 
from the wave-vectors. 

Let us consider the simplest situation with the three vectors Na aligned along a given 
direction of space with spherical coordinates {9,4)). It is easy to show that in this 
situation all of the coefficients lap^y are equal to zero, so the bispectrum contribution 
from (50) becomes zero as well. Notice that this does not apply either to the power 
spectrum or to the Abelian part of the bispectrum, which retain a non-zero anisotropic 
part (see Eqs.(15) and (30)) when the vectors Na are parallel to each other. 
There are several other situations that can be taken into account, for example one of the 
Na being aligned along or perpendicular to one of the ki and the other two generically 
oriented, two of the A^^^ aligned along ki and k^ respectively and the third one generically 
oriented and so on. Our next step is to provide a plot of the bispectrum complete of its 
isotropic and anisotropic parts in the favoured configuration (in our case the local one). 
For convenience while plotting, let us normalize the vector Na so that they all have 
the same length and let us restrict our attention to situations where the total number 
of angles the la/sj coefficients depend on does not exceed two. One possibility is for 
example given by the following case 

iV3 = iV^(0,0,l) (70) 
Ni = N2 = A^A(sin 6* COS0, sin 6' sin 0, cos6'), (71) 

where, as mentioned above, we are working in the coordinate frame with k^ = x and 
ki = k2 = z and we took Na = \N"-\ for all a = 1,2,3. Let us introduce the angle S 
between A^i 2 and k^. Then from Eq.(50) we have 

Bq D gl—^ ^"/37(^' 5)Rci3^{x\x2, X3). (72) 
^1 a/37 

The functions Raf^-yix^, Xs) are plotted in Fig.l for x* = 1 and their analytic expressions 
are provided in Appendix C. See Appendix D for the expressions of the coefficients 

A plot of the 'non- Abelian' bispectrum normalized to the ratio {g^H'^rn^ Nj^ / [k^x^x^) 
is given in Fig.2 for fixed values of x*, X2 and x^ (see Appendix D for its analytic 
expression) . 

6. Calculation of the non-Gaussianity parameter f^^ 

The non-Gaussianity of a given theory of inflation and cosmological perturbations can be 
studied by looking at the expression of the non-hnearity f^L- This parameter is deflned 
by the ratio of the bispectrum to the square of the isotropic part (as for example given 
in Eq.(18)) of the power spectrum 
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Figure 2. Plot of ,5) = <5, x*, 0:2, X3)a;ia;ifc6)/(g2^2^,2^4)]^^^^^^^^^^^ ^^^^ 

As wc can sec, the graph dependence from 9 is extremely weak; in fact, the dominant 
contribution from Eq.(D.9) can be written as f{9, S) = —140 cos^ S. 



5 Pi'°{ki)Pi'°{k2)+ perms.' ^ ' 

As we saw in the previous sections, the complete expression for the bispectrum is 
made up of an 'Abehan' and a 'non-Abehan' parts, which can be read respectively from 
Eq.(21) and (50). We are going to rewrite their sum below in a compact expression 



B^{kr, h, h) = Bfih, h, h) + ^^(^1, h, h)- (74) 

Our next step will be to calculate the contributions from the two terms on the right hand 
side of the previous equation to the f^L parameter and establish a comparison between 
them. The 'Abelian' part of the bispectrum was given in Eq.(21) as the contribution of 
three different pieces 



B^{ku ^2, k,) = B^^ + i?f + B^\ (75) 



where 



Bf' = {N^f N^^P\ (76) 

B^^^{NAfNAAP\ (77) 

5f ^ N^NaN^aP\ (78) 

where the symbol = stresses the fact that we are interested in the overall order of 
magnitude of the quantities above in terms of the parameters of the theory, without 
worrying about the details of numerical factors of order one. Also, we now indicate the 
power spectrum of the scalar field perturbation by P. At late times this is proportional 
to the power spectrum of the gauge fields. This can be seen by looking at the transverse 
mode function for 5B in Eq.(57) and at the longitudinal mode in Eq.(58), which, after 
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taking the limit x — » 0, give P+ = P and Piong = 2P. We shortened A^* by A^^, A''^-^ by 
Naa and A^^^ by A^^^ in order to make the expressions simpler, suppressing the gauge 
and the vector indices. Notice that in the theory we have chosen, no interactions occur 
between the scalar and the gauge fields, therefore we can set A'',^^ — 0. 
Similarly, we rewrite the 'non-Abelian' part of the bispectrum in a simplified fashion 
that leaves out numerical factors and complicated functions of the external momenta 

B^^{k„h,h) = 9'cH'A{NAf, (79) 

where we now call A the background value of the spatial part of the gauge fields 
(suppressing all indices). 

Let us now rewrite the isotropic part of the power spectrum 

p-°^ (A^^)2p(i+/5), (80) 
where f3 = {NA/N^f. 

Before proceeding to the evaluation of f^L, we need to make more assumptions about 
the scenario we are dealing with, since the derivatives of A^ are model-dependent. Two 
possible choices are the vector curvaton model (introduced for the first time in [31], see 
also [34] for one of the possible realizations of this paradigm) and vector infiation [33]. 
The former is the analogue of the (scalar) curvaton mechanism: a vector field plays the 
role of the curvaton, i.e. a field whose energy density is subdominant during infiation 
and, after the infiaton has decayed, it is is mainly responsible for the production of 
the curvature perturbation. In vector infiation, the universe exponential expansion is 
driven by vector fields, which need to be spatially oriented in such a way as to avoid 
large scale anisotropies; one possible field configuration is a triplet of orthogonal vector 
fields, another one consists in a large number of randomly oriented vector fields. In our 
model we are dealing with a triplet of vector fields, so if we assumed the fields were 
orthogonally oriented w.r.t one another, vector infiation would be a possible scenario 
to work in. A third possibility can be to consider a scalar field driven infiation in the 
presence of the SU (2) gauge multiplet and to assume that the latter is undergoing slow- 
roll and its energy density is subdominant w.r.t. the energy density of the scalar infiaton. 
This way we keep the vector fields in the game, without imposing any restriction on 
their orientation in space and, at the same time, we avoid unwanted anisotropies in the 
power spectrum. We are going to consider both this last possibility, which we dare refer 
to as 'vector infiation' for simplicity, and the curvaton scenario. In 'vector infiation' A"* 
is given by Eq.(B.12), from which we get 



In the curvaton scenario [17] 
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where r = (3pa) / (Sp^ + 4p(^) , and being respectively the energy density of the 

vector fields and of the infiaton at the time where the vector field(s) decay. 

We are now ready to evaluate /jv^ as the sum of three contributions corresponding to 

75 tree 75 tree „„J 75 tree 

Jnl = Inl + Inl + Inl^ (85) 

where 

1 AT 

(86) 
(87) 
(88) 

Numerical coefficients of order one were not reported in the previous equations, m is a 
quantity with the dimensions of a mass and proportional to the Planck mass in the case 



fM _ 
Jnl — 


1 N^^ 


(1 + m ' 


fM _ 
Jnl — 


/3 Naa 




pNA 

Jnl — 


f3^ 2 fm 


(l + /3)2^^ [h 



of vector infiation and to Atotj (where we define Atot = \/J2b from Eq.(84)) in 
the vector curvaton model. 

We summarize our results for the two models in Table 1, where we used the expressions 
2^ (mpy/e^)""^ and N^^ 2^ nT'p^, with e^j, = {(p"^) / {2mpH^) . The background values A 
and Atot together with and H are all meant at time x = x*. 



Table 1. Order of magnitude of f^L in different scenarios. 



Jnl 



fA2 

J NL 



fNA 
J NL 



.infiation 



mpH 



V. curvaton 



9 



Am?p 

^tot 



/ Amp \ ' 



A^m% 



The ratios between the non-Abelian and the Abelian parts of the /tvl parameter 



are 



NA 

fAi — *'<t>yc 

J NL 



A 



H) \mp ) 



fNA . A 

fA2 — yc \ TT 

J NL 



H 



(89) 
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in vector inflation and 



re^ \h) ' f^l r^e^ \H ) \mp) 



for the vector curvaton. 

When estimating the order of magnitude of /tvl as calculated above, it is important to 
remember that our results in Eqs. (86) through (90) and in Table 1 have been derived 
neglecting all vector and gauge indices. Also, these expressions involve the background 
value of the slowly- rolling flelds A'^ in a multifleld space and so f^i depends on the 
whole speciflc background conflguration. In this respect, different possibilities can be 
considered and it turns out that, depending on the chosen conflguration, the non-Abelian 
contributions to /tvl can subdominant or dominant w.r.t. the Abelian ones. Notice 
that the ratios QcA/H = {gcA/mo) [niQ/H) for the different gauge flelds are bounded 
from above, as shown in Eqs. (A. 9) through (A. 11). There is a very large region of 
the parameter space of the background gauge fleld conflgurations for which this upper 
bound is much larger than one and, therefore, it does not prevent f^^ from dominating 
over f^i^. Here is an example: suppose that two of the gauge flelds, for instance and 
B are aligned, about equal in mag nitude {\B^\ ~ l^^l) and that \B^'^\ > \B'^\. In this 
case, by looking at Eqs. (A. 9) and (A. 10), it is easy to realize that the upper bounds 
can be much larger than one. Therefore, the ratios in Eqs. (89) and (90) can be much 
bigger than one depending, for example, on how A/H compares to the other quantities 
appearing in these expressions (notice: the theory puts no constraints on A/H which 
can in principle be a very large number). Obviously, inspecting Eqs. (A. 9) through 
(A. 11), if the A"'s are not approximately aligned and their magnitude A is roughly the 
same, then the upper bound gc <^ H/A holds. 

As to the absolute value of /jvl > expanding around /? -C 1, from Eq.(88) we have 

respectively for vector inflation and the vector curvaton model. The ratio mp/H is of 
order 10^ and A/H could be much bigger than one as well, so in principle f^^ can be 
much larger than one in both models. Equivalently, from Table 1, f^^ can be put in 
the following form (as a function of the slow-roll parameter e^) 

fNA _ I 9^4 (^)^ (^1)' ^ vector inflation 

g^elr^ (^) {^) vector curvaton 

7. Overview and conclusions 

We have considered a triplet of SU (2) vector bosons non-minimally coupled to gravity, 
Eq. (4), in order to have ~ —2H^. Using the 5N and the Schwinger-Keldysh for- 
malisms we have computed the contributions to the curvature perturbation three-point 
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correlation function arising from the gauge fields self-interactions. These interactions are 
of two kinds, third-order interactions (proportional to one power of the SU{2) coupling 
gc) and fourth-order ones (proportional to g^). The former provide a vanishing con- 
tribution to the bispectrum because of the antisymmetric properties of the Levi-Civita 
tensor appearing in the Lagrangian, the result being independent of the particular form 
of the wave functions appearing in the gauge bosons operator expansions. The quartic 
interactions produce instead a non-zero contribution, Eq.(50). This result is anisotropic 
if the wavefunctions of the transverse and longitudinal modes of the vector fields are 
different from each other, which seems to be the case by looking at their equations of 
motion, Eqs. (A. 16) and (A. 17). 

The ongoing debate about the instability of the longitudinal mode sheds some 
doubts about the true physical meaning of the mode. Our analysis, in this respect, is 
pretty general, being mainly focused on studying the non-Gaussianity effects in non- 
Abelian theories. In addition to that, we tried to also use a more general approach 
when dealing with longitudinal modes, with a parametrization of the wavefunctions in 
terms of the transverse mode and of an unknown function n{x), Eq.(68), which, to a 
first approximation, can be set equal to a constant near horizon crossing, n{x) = n{x*). 
A particular case of this parametrization occurs when the longitudinal and transverse 
modefunctions coincide, which offers an interesting result, i.e. the isotropization of the 
full bispectrum, both the 'non-Abelian' and the 'Abelian' parts, Eq. (69) and (21), re- 
spectively. We are aware that the particular Lagrangian we used represents only one 
possible model of primordial vector fields. Other theories have been proposed, which do 
not present any kind of instability (see for example, [55, 36, 35, 56]). 

In the case of an anisotropic bispectrum we showed that it can be written in terms 
of isotropic parts, i.e. functions of the moduli of the external wave vectors ki, modu- 
lated by anisotropic coefficients la/s-y which depend on the angles between the (wave and 
gauge) vectors. We studied the shape of the 'isotropic parts' of the bispectrum, Eq.(50) 
and Eq.(59) through (62), which turned out to peak in the local momenta configuration 
(see plots in Fig.l). Using this finding, we analysed the 'anisotropic' part, Eq. (63) 
through (66) for different spacial configurations of the gauge vectors. A limiting case is 
represented by the three components of the SU{2) gauge group being all aligned with 
one another: because of the presence of the Levi-Civita tensors and with simple sym- 
metry consideration, the contribution from (50) to the bispectrum in this case is proven 
to be zero. This result could be interpreted as the Abelian limit of our theory, since 
the three components of the SU{2) multiplet are all identifiable with a unique spatial 
direction. 

Another example of vector fields configuration was provided in Eqs.(70)-(71). For this 
particular complete plot of the bispectrum was given in Fig. 2 for the local con- 

figuration. It contains the information about the anisotropy of the bispectrum, showing 
how its amplitude is modulated according to the angles between the wavevectors ki and 
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the preferred directions induced by the vector fields. 

Finally we have calculated the parameter Jnl, estimating the level of non- 
Gaussianity. We analyzed the different contributions to /atl, i.e. the 'Abelian', Eqs. (86) 
and (87), and the 'non- Abelian' ones, Eq.(88). We compared them both for vector in- 
fiation and for the vector curvaton models and for different background configurations. 
It turns out that one contribution can be dominant w.r.t. the other or viceversa, de- 
pending on the given background configuration, Eqs. (89) and (90). Focusing on the 
order of magnitude of the 'non-Abelian' contribution, we noticed that it can be much 
bigger than one, for a large region of parameter space, Eq.(91). 

In a forthcoming paper [54], we will study the trispectrum of the curvature 
perturbation in the same vector fields populated model, calculate its magnitude and 
shapes and and investigate the relationship between the f^i and tnl parameters of the 
theory. 
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Appendix A. Background and first order perturbation equations for the 
gauge fields 

The equations of motion for the gauge fields have been completely derived for the U{1) 
case in [31]. We are going to carry out a similar calculation for the SU{2) case 



2 
I' 

The = component of the equations of motion is 



(A.l) 



where Fj;^^)"^ = d^B^ - d^B""^ 



+g,e'''''B'^B'^B^'] = 
where B^ = B^{x,t). 

Let us now move to the spatial (z/ = i) part of (A.l) 



d,B']B?, - 2B''diB^ - B'^B^ 



3 3 



B1 + RBI - \djdjBl + M'Bl - d.B^ - RdiB^ + \d.^,B] 



(A.2) 



abc 



HBlBt 



+9ce 



abc 



diBf^WQ — !?■ B, 



b rtc 



i -^0 



-abc 



{d,B])Bt + B]d,B\ 
(d.B^B^-(d,B^Bi 



, 2 abc bb'c' 

-j-g^b c 



rjft' rjC 



_ 9c abc bb'c' 



B'^B^'Bf 



(A.3) 



If we contract Eq.(A.l) with d,^, we get the integrability condition 

{aMYB^ - M^diBf + 3H[dABS - d^Bf) + g,e'''''[2H[d^B^B'Q + B'^B'o + B]Bp^ - d^B'^B] 
-B)B'^ + d.BlB'^ - d^B'oB'o - diB^B'^ + ^(B^d'^B'r + diBp^B^ + 5^9,9^5,'= + d^B]B''^ 



-d,d,B\B^^ 



2 abc bb c 
9 c ^ 



I IDC ryb r)C , rye r>b rye , rye jyb jyc 



c Tyb 



0-^0 -^0 



2iJfi,^5o^ Bt 



RC Tyb jyc uc jyb rye f\ jyb jyc JDCf\ jyb rye jyb _n jyc 

= 



+-[diB^BfB^' + Bp.BfB^' + Bp^'diB;: 

which reduces to Eq.(7) of [31] in the Abelian case. 
Combining Eq.(A.4) with Eq.(A.2) we get 



(A.4) 



{aMfBg - W'diB^ + 3H[a'M''B^ + gce'""[ - - '^B'jOjB'q - IB^B'^ + g.e"'''' Bpl B 

+g,e''''[2H(diBtBl + B^B^ + Bpf) - djB'.B^ 

-i3)B'^ + d.BlB'^ - d^B'^B^ - d^B\B'^ + ^iB^d^B^ + diBpjB^ + Bp.diB'r + d^B^B''^ 



?b' Dc' 

j 
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I 2 abc bb c 
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d,B^,B\ B^, 



B^diBi Bq - BqB^ diB^ 



+-[diBpfB^' + Bp.B^B^' + BpfdiB^' 



(A.5) 



Plugging this into Eq.(A.3) we get 
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(A.6) 



Let us consider the background part of the vector fields, i.e. diB'^ = 0. Then from 
Eq.(A.2) 



a^M^B^ + g.e'''" 



B^B] + g,e'''''''B]B'^Bf 



0. 



(A.7) 



Before proceeding with the derivation of the equations of motion for the background and 
the field perturbations, it is necessary to make some comments on Eqs. (A.6) and (A.7). 
One approximation that we have been using in this paper is allowing the A" fields to 
undergo slow-roll during infiation. One possible way to achieve this is by restricting the 
parameter space of the background gauge fields through the request that their temporal 
components should be much smaller than the spatial ones, Bq <^ \B\\/a{t), h = 1,2,3, 
and, in addition to that, assuming B^ ~ Bq, b,c = 1, 2, 3. With these assumptions, 
the temporal component can be factored out in Eq. (A.7), using the approximation 
B^ ^ HB1 (valid in a slow-roll regime). A solution to (A.7) is then given by Bq = 0. 
Adopting this solution and plugging it in Eq (A.6), it is easy to show that a slow-roll 
equation of motion for the physical fields 



follows from (A.6) if H <^ rriQ and 







(A.8) 
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HA'] 


^ - {A^f cos2 


^'23 


- (A1)^C0S2 



(A.9) 
(A.IO) 
(A.ll) 



are satisfied. In the equations above, we defined A"- = \A°'\ and cosOab = A"- ■ A'', a 
and b running over the gauge indices. The quantities appearing on the right-hand sides 
of Eqs.(A.9) through (A.ll) can be either large or small w.r.t. one, depending on the 
specific background configuration, i.e. on the moduli of the gauge fields and the angles 

Suppose now the conditions described above are all met, then from Eq (A. 6), in terms 
of the comoving fields, we have 

Bf + HBf + M'B^ = 0. (A. 12) 

Let us now derive the equations for the perturbations. Eq.(A.2) becomes 



dj6B] - d'6B^ + a'M'6B^ + g^e 



abc 



d,5B]B'^ - 2Bp,5B'^ - 5B'jB^ - BpB^ 



I I 

cb c 



5B]Bl B^ 



B]5Bl fij 



(A.13) 



Eq.(A.l) for the field perturbations gives 



55f + H5i3'^ - \djdj5Bf + M^SB^ + \didj5B^ - HdiSB^ - HdiSB^ 



9c_ abc 



d,5B])Bl 



Bp,5Bl 



9c abc 



d,bB'^B]-[d,bB\)B] 



_9c_^abc^bb'c' 



+9ce 
+9c£ 



abc 



dBp^ 5; + B^6B', + Bp'i 6B^ 



H [B^dBt + 5B',Bt) + 5BlBt + B^dB^ + 5B',Bi 

B'^SB' 



abc 



d,5Bl5Bl - 6B^B^ 



: 2 abc bb c 
'9c£ £ 



5BIB\ 5S 



B15B\ B^, 



B^,B\ 5BI 







Finally from Eq.(A.6) we get 



Bl5Bt 



(A. 14) 



5Bf + H5Bf - \d'SB^ + M^SBf + 2Hd,6B^ + (~ g^terms) = 0. (A. 15) 

When calculating n — point functions for the gauge bosons, the eigenfunctions we need 
are provided by free- field solutions, i.e. by solutions of Eq.(A.15) with g^ being set to 
zero. This is exactly the Abelian limit, in fact in this case Eq.(A.15) corresponds to 
(18) of [31] and can be decomposed into a transverse and a longitudinal part 
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(9^ 



2k^ 



SB' = 



k^ + {aMyj 

where the time derivatives are intended w.r.t. cosmic time. 



= 



(A.16) 
(A.17) 



Appendix B. Calculation of the number of e-foldings of single- (scalar) field 
driven infiation in the presence of a vector multiplet 

Let us consider the complete Lagrangian of our theory as in Eq.(4). Let us assume that 
the SU (2) gauge multiplet undergoes slow-roll as well as the scalar field but the latter 
provides the dominant part of the energy density of the universe. This last hypothesis 
is necessary in order to produce isotropic infiation (i.e. in order for the anisotropy 
in the expansion that the vector fields introduce to be negligible w.r.t. the isotropic 
contribution from the scalar field) . The expression of the number of e-foldings is 



N = N scalar + 



vector 



^scalar ~(" 



4"^P a=l,2,3 



J2 A^-A^ 



(B.l) 



The previous expression can be easily derived from the equations of motion of the system 
neglecting terms that are proportional to the SU (2) coupling constant gc and assuming 
slow-roll conditions for both the scalar the gauge fields. 
The starting point is represented by Einstein equations 



scalar ~l~ Pvector^ 



(B.2) 



where we split the energy density into a scalar and a vector contribution. In slow-roll 
approximation, p scalar ~ ^(</')- Let us calculate Pvector- The energy momentum tensor 
for the gauge bosons 



T. 



vector 



6L 



5g 



(B.3) 



where, as a remainder, L = -{l/A)g^"^g''f^ F^^F^p + {M"^ /2)gf"'B1^B^. So we get 



rpvector 
-^00 



2a2 



"^0 na-na 
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^BtBt 
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2a2 



^abc^ab c J^^ J^'r Q'r 
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BfBf 



(B.4) 



where sums are taken over all repeated indices. Let us write this in terms of the physical 
fields 



2 2 2 

2 

+ ^e'^^'^e'^^'^'AM^f A?' (B.5) 

If we neglect the non-Abelian contribution and we set = 0, we are left with the 
Abelian result [33] 

.^.ector = ^ + -O^M" (B.6) 

— * 

The equation of motion for the background vector multiplet A"- can be derived from 
Eq.(A.12) 

+ 3i/i^ + mlA^^ = 0. (B.7) 

which is equal to the equation of a light scalar field of mass mo, if mo <^ H. If the 
conditions for accelerated expansions are met, Eq.(B.7) reduces to 

+ mlA'l ~ 0. (B.8) 

We are now ready to derive Eq.(B.l). Let us start from the definition of N and keep 
in mind Eq.(B.2), where we are assuming the existence of a scalar fields (f) in de-Sitter 
with a separable potential governed by the usual (background) equation 



+ 3i70 + 1/'=O (B.9) 
and slowly rolling down their potential. Then we have 

N = /' Hit = /' H-'^i = 8.G /' ^-Mii + /' ^-^,lt 
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where A'' = A" ■ A''. So 
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hit*) 3H(f) aJA<^it') V 4 / 3HA';A1 

1 r<l>{t) V 1 /m2\ AM? 

= - — TT^^*^ + — E / \-r 
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after using the slow-roll conditions. Eq.(B.l) is thus recovered. 
In the final expression for the bispectrum then we can substitute 



K^^=(^]a^ (B.12) 

where the derivatives are as usual calculated at the initial time if . 
The upper limits in integrals such as the ones in Eq.(B.ll) depend on the chosen path in 
field space and so they also depend on the initial field configuration. It is important to 
notice though, as also stated in [53] , that if the final time is chosen to be approaching (or 
later than) the end of infiation, the fields are supposed to have reached their equilibrium 
values and so N becomes independent of the field values at the final time t. Eq.(B.12) 
is thus recovered. 



Appendix C. Complete expressions for the functions appearing in the 
bispectrum from quartic interactions 

We list below the complete expressions for the functions -Ba/37 and Ca^^ (a, /5, 7 = 

R,V) appearing in Eq.(59) through (62) 
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where k = ki + k2 + k^. 



Appendix D. Profile of tlie bispectrum: detailed expressions of the plotted 
functions 

We are going to study the profile of the 'non-Abelian' contribution to the bispectrum 
in terms of the angles between gauge and wave vectors in the local configuration and 
for the set up considered in Eq.(70) and (71). The coefficients lap'y then become 
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where m? = (A)/{Na), A being the background value of the A^s evaluated at horizon 
crossing. 

Fig. 2 shows a plot of the 'non-Abelian' bispectrum normalized to the ratio 
{g'^H'^m'^N^)/{klx2xl), as a function of the angles 6 and 6 and for fixed values of x*, X2 
and X3 



B^{e,S) r:^gl—^——^ cos^5(8cos^- 1.4 X 10^) + 3 cos5(cos'^ ^ - cos^ ^ - 11) 
- 11 cos 25 - 40 - 6 cos 20 - cos ^(3 cos^ ^ - 30 cos ^ - 10)" 
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[D.9) 



where we set x* = 1, while X2 and X3 were chosen in the 'squeezed' region, X2 = 0.9 and 
X3 = 0.1. 
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